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Abstract 



^ ' Number theorists have studied extensively the connections between the 

JiJ . distribution of zeros of the Riemann (^-function, and of some generahza- 

(^ ', tions, with the statistics of the eigenvalues of large random matrices. It 

^^ I is interesting to compare the average moments of these functions in an in- 

Q>^ ' terval to their counterpart in random matrices, which are the expectation 

0^ ' values of the characteristic polynomials of the matrix. It turns out that 

,^ . these expectation values are quite interesting. For instance, the moments of 

order 2K scale, for unitary invariant ensembles, as the density of eigenvalues 
j^ I raised to the power K^ ; the prefactor turns out to be a universal number, 

J3 ' i-e. it is independent of the specific probability distribution. An equivalent 

behaviour and prefactor had been found, as a conjecture, within number 
^ . theory. The moments of the characteristic determinants of random matrices 

k> , are computed here as limits, at coinciding points, of multi-point correlators 

H ' of determinants. These correlators are in fact universal in Dyson's scaling 

limit in which the difference between the points goes to zero, the size of the 
matrix goes to infinity, and their product remains finite. 
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1 Introduction 

The correlation function of the eigenvalues of large N x N matrices are known to 
exhibit a number of universal features in the large- iV limit. For instance in the Dyson 
limit |I], 0], when the distances between these eigenvalues, measured in units of the 
local spacing, becomes of order 1/A^, the correlation functions, as well as the level 
spacing distribution, become universal, i.e. independent of the specific probability 
measure. For finite differences, upon a smoothing of the distribution, the two-point 
correlation function is again universal f^, Q. The short distance universality was also 
shown to extend to external source problems 0, ^, 0, |^, in which an external matrix 
is coupled to the random matrix. 

In this article, we study the average of the characteristic polynomials, whose 
zeros are the eigenvalues of the random matrix. The probability distribution of the 
characteristic polynomial det(A — X) of a random matrix X, a polynomial of degree 
A^ in A, may be characterized by its moments < det (A — X) >, or better by its 

K 

correlation functions < JJ^ det(Ai — X) >. 

1=1 
This study is motivated by various conjectures which appeared recently in number 

theory for the zeros of the Riemann (^-function and its generalizations known as L- 

functions |]T^. Indeed the characteristic polynomials, as well as the zeta-fuctions, have 



their zeros on a straight line, and these zeros obey the same statistical distribution. 
For the 2i^-th moment of the Riemann ^-function {K is a positive integer), it has 
been conjectured 0, |T^ that 

^ r dt\a\ + tt)\''' ^ 7xa,.(logT)^' (1) 

I Jo 2 

where ax is a number related to the Dirichlet coefficient (the divisor function) dxin), 

and 

The explicit formula for ax is given in the Appendix, together with summation 
formulae for the Dirichlet coefficients, which are related to (|1|). In this work we shall 



compute the equivalent of (Q) for random matrices, show that the density of states 
p(A) replaces logT, and that the same number 77^ is universally present. 

For the negative moments, similar conjectures have been proposed, with a cut-off 



parameter 6 for avoiding divergences |Tl[], and we show here how to obtain these 
negative moments for random matrices. 



Several types of L-functions have been introduced [|T2|, which correspond to the 
three standard classes of random matrices The conjecture for the average of the 
moments (|l]) has been extended to these L-function [Q. The average is taken as 
a sum of the discriminant d, for instance, for the Dirichlet L{^,Xd) function. The 
relations between the distributions of the eigenvalues of the random matrix theory 
and the statistical distribution of the zeros of the various L-functions has also been 



studied p|, |T|]. 



Our aim in this article, is to clarify the universality of the moments of the charac- 
teristic polynomials for these three classes. The circular unitary ensemble, has been 
studied earlier by Keating and Snaith[|10|, who did obtain the 7^ in (0) from their 



calculation. However this ensemble has a constant density of states, and furthermore 
it does not allow to study the universality of these properties. In this work we have 
considered a Gaussian ensemble and non-Gaussian extensions, instead of the circular 
ensemble, to verify both the explicit dependence in the density of states and the 
universality of the coefficient 'Jk- In the process of the derivation, we have found it 

K 

necessary to start with the K-point functions < ]^det(A/ — X) >, which are shown 

1=1 
to be themselves universal in the large-N Dyson limit, in which N{Xi — Xj) is held 

fixed. The moments are then simply the limit of these functions when all the Dyson 

variables vanish. 



2 Correlation functions of characteristic polyno- 
mials 

We consider random M x M Hermitian matrices X with a normalized probability 
distribution 

P(X) = |exp-iVTrV(X), (3) 

in which ^ is a given polynomial. It will turn out to be convenient to distinguish 
here M and A^, but we later restrict ourselves to a large N and large M limit, with 
liuiM/N = 1. Let us consider the correlation function of K distinct characteristic 
polynomials : 

K 

Fx(Ai,---,Ax) = (ndet(A,-X)), (4) 

in which the bracket denotes an expectation value with the weight (0). 
Integrating as usual over the unitary group, we obtain 

1 „ Af K M 

fk{Xi, ■■■,xk) = ^^ / n dK^i) ^^{^i,---, ^m) n ik-^" ~ ^i) (5) 

in which dii{x) denotes the measure (i/x(x) = dxexp—NV^x), A the Vandermonde 
determinant A(xi, ■ ■ ■ , xm) = H (^^i ~ ^j)> ^^^ ^m the normalization constant 

l<i<j<M 

„ M 

Zm = Yl dfJ'i^i) ^^(a^i, ■ ■ ■ , Xm)- (6) 

•^ 1 

We now use the obvious identity 

A(a;i,---,XM) 11 IKAa-Xi) = — , (7) 

a=ii=i I\{Xi,---,\k) 

and represent the Vandermonde determinants A(xi, ■ • • , xm) and A(xi, • • • , xm; Ai, • • • , Xk) 
as determinants of arbitrary polynomials whose coefficients of highest degree are equal 
to unity (the so-called monic polynomials) 

Pn{.x) = x^ + lowerdegree. (8) 



Then 



A{xi,---,xm) =detpnixm) 
[n runs from zero to M — 1 and m from one to M),and 

A(xi, ■ ■ ■ , xm; Ai, ■ ■ ■ , \k) = detpaiub) 



(9) 



(10) 



in which a runs from zero to M + K — 1, b from one to M + K and Ub stands for Xb 
iib<M,oTXbioTM<b<M + K. 

Choosing now the polynomials orthogonals with respect to the measure dfj, : 



:ii) 



Pn{x)pm{x)dfi{x) = hnSnr 

we may easily integrate over the M eigenvalues 

» M /A/-1 \ 

/ Yldjj,{xi) A{xi, ■ ■ ■ , xm; >^i, ■ ■ ■ , >^k)M^i, ■ ■ ■ , xm) = M\ i Yl hnj det pa{X^), 

(12) 
in which a runs from M to M + i^ — 1 and jS from 1 to K. Similarly the normalization 
factor Zm is given by 

M /M-1 



„ M /M-1 \ 

Zm = J liM^i) ^^^u- ■ ■ ,xm) = Ml i UhA. 



(13) 



We thus end up with 



Fk{Xi, ■ ■ ■ , a 



1 



K 



A(Ai,---,A 



■det 



K) 



Pm(Ai) PAf+i(Ai) 

PAf(A2) PM+i{X2) 



Pm+a'-i(Ai) 

PM+K-i{X2) 



.(14) 



Pm{Xk) Pm+i{Xk) ■■■ Pm+k-i{Xk) 
If we are concerned simply with the moments of the distribution of a single charac- 
teristic polynomial, we obtain from ([T^) 



/Xi^(A) 



F^(A,---,A) = ([det(A-X)]^) 

PmW Pm+i(A) 
p'mW p'm+iW 



(_l)i^(^-i)/2 



det 



PAf+i4'-l(A) 

Pm+a'-i(A) 



.(^-1)/ 



.(^-1) 



(K-l) 



p'm'W P)vtT(A) ■■■ p'm^kUX) 



(15) 



These expressions are all exact, but in the next section we shall be concerned with 
the large N limit. Then (i) the interesting case is that of even K, since for odd K 
the result is oscillatory ( for instance ioi K = 1 fii{\) = PmW ), (h) it will turn out 
that, even if we are interested simply in the moments fixW, it is more convenient 
to study first the large iV-limit of the F^ with distincts Aj and afterwards let them 
approach a single A. 

The results that will be derived later for those -Fr-'s and /ii^'s will be shown to 
be universal in the Dyson limit, in which A^ goes to infinity, the Aj — Aj goes to zero 
for any pair i,j, and the products A^(Aj — Xj) remain finite. We first derive explicit 
formulae for the Gaussian case, and show later that they do apply to any random 
matrix distribution P{X) of the form @. 

3 The Gaussian case 

We now specialize the result (|T^) of the previous section to the Gaussian distribution 
of M X M Hermitian matrices 

P(X) = ^exp-^TrX2, (16) 

with 

M = N-K, (17) 

(the reason for this choice of M will be clarified in the next section). Then the 
polynomials that we have introduced, are Hermite polynomials, and with our nor- 
malizations, 

Hn{x) = ^e^^^/^(^)"e-^^^/2 = X" + l.d., (18) 



and 



K = — \\ — - 19 



The integral representation 



5 



over a contour which circles around the origin in the z-plane, turns out to be well 
adapted. A repeated use of this formula in the result ( |T4|) yields 



F2k{^1, ■ ■ ■ , ^2k) 



2/r / ^ \ 2K 2 



We can expand the determinant in the r.h.s. and keep only one of the {2K)\ terms, 
antisymmetrizing instead the integration variables zi. This gives 



'"^ ^ dzi \ , ^^^,zf 



This expression for the expectation value of a product of 2K characteristic polyno- 
mials, as an integral over 2K complex variables, is exact for finite A^ and M. 

We are now in position to study the large iV-limit through a saddle point in- 
tegration over each zi. Since we have chosen M -|- A' = A^ each z has a weight 
— rr exp —N{z'^/2 + Xz + log 2;), which presents two saddle points z±, solutions of the 
equation z"^ + Xz + 1 = 0, i.e. with the parametrization 

A = 2sin0, (23) 

when A lies on the support of the asymptotic Wigner semi-circle of the density of 
levels, 

z+ = ie'^, z. = -ie-'^. (24) 

Therefore there are, a priori 2^^ saddle-points at which the moduli of the weight 

2K ^2 2K ^2 

exp [— A^y^(— + XiZi + log -2;)] are the same. However, it is only when ^^{— + 

1 2 I 2 

XiZi + \og Zi) is real (in the Dyson limit in which the differences between the A's are 
small), that the oscillations, which damp the result, are not present. Therefore we 



keep only the j^ saddle-points in which we take K solutions of type 2;+ and K 
of type Z-. We are now interested in Dyson's short-distance limit. Defining 

1 2K 

A = ^EA., (25) 

and the density of eigenvalues at this point 

p(A) = -lV4^A2 = -cos0, (26) 

/vr TT 

we introduced the scaling variables 

2K 

X, = 27rArp(A)(A„-A),with5]x„ = 0, (27) 

a=l 

which are kept finite in this limit. Then the fluctuations around a saddle-point may 
be taken all at the point A, and they yield a factor 

(|)^[(l-4)(l-.!)]-^/^ = (iVp(A))-^ (28) 



We must now take into account the various factors in ( p2D at these saddle-points. 

2K 

In the Dyson limit the factor J^ z^ which remained in the denominator, may be 

1 
replaced by one, since at a given A one has z+z_ = 1. The only delicate factor is 

Afz ■ ■ ■ z k) "^^ z"^ 

thus ' ^ — rexp [—N^^{— + A/Z/ + log^^)], which we must first compute at 

A(Ai, ■ ■ ■ , X2k) 1 2 

one of the saddle-points, and then take the sum over the ^ j saddle-points. We 

consider first the saddle-point 

zi{Xi) = z+{\{) / = 1,---,K 

zi{Xi) = z^{\i) l = K + l,---,2K. (29) 

2K ^2 

If we expand in the Dyson limit the weight exp [■~N'^{ h XiZi + log^;^)] one finds 

1 2 

■2K -1 

. Z-] 



exp[-ArE(^ + Ai2;z + log2;0] 
1 ^ 



A^ 
exp A^K(1 + ^) X exp -A^ 



K 2K 

5:(A,-A)^+(A)+E(A.-A)^-(A) 

1 K+l 



, (30) 



(we have used — (-z|(A) + \z± + \ogz±) = z±). Therefore at that saddle-point, in 

UA 2 



terms of the scahng variables (27) 



2K 2 \2 K 

exp [-NJ2{^ + Xizi + log zi)] = expNK{l + — ) exp -tj^xi. (31) 

Let us consider now the ratio of Vandermonde determinants at that same saddle- 



point: 



A(zi,---,2:2i^) _ -j-j- Z+(A;) -2:+(Am) -r-r Z-(A;) -2:_(An,) 

A(Ai, ■ ■ • , \2k) X<l<m<K ^l " ^™ K+l<l<m<2K ^' " ^"^ 

Z+{Xl) - Z^{Xm) 



n "^ , _r^ '"^ (32) 



In the scaling limit, this factor becomes 

A(zi, ■ ■ ■ , Z2k) ,dz+dz. 



^(^-i)/2(2zcos</))^' n 



A(Ai,---,A2ii-) dX dX i<i<K,K+i<m<2K -^i ~ -^m 

= {mf\2^p{X)r^^' n -4:-- (33) 

l<l<K,K+l<ni<2K "''' ^»^ 

Leaving aside for the moment the overal factors which do not change at the various 

saddle-points, we note the result from this particular one which is 

exp —i y^ xi Y\_ , and consider summing over the ^ ] saddle- 

1 l<l<K,K+l<m<2K ^l ~ ^m \ J 

points. The sum is best done under the form of an integral over K variables. Indeed, 
if we consider 

over a contour in which each Ua circles around the x's, we recover exactly the contri- 
bution previous saddle-point by choosing ui = xi, ■ ■ ■ ,uk = xk, or any permutation 
of those K x's. In view of the Vandermonde in the numerator, all the u's have to 
be different, and thus there are indeed ^ ] poles to be added, which reconstruct 
exactly the sum on the saddle-points that we needed to perform.. 

Collecting the various factors that came on the way, we end up with the final 
formula 

exp-{-J2V{Xi))F,K{Xi,---,X2K) = 
^ 1=1 



(2.NpiX)) ^^ fn^-^p -(E «.) ng..n?A(.„-x,) - P'' 

If we specialize to i^ = 1 one finds 

/v Sin T" 
exp{--(\/(Ai) + \/(A2))}F2(Ai, A2) = e"^ (27riVp(A)) (36) 

^ Jb 

witli X = 7rA^p(A)(Ai — A2), in which we recover the well-known Dyson kernel, which 
characterizes the correlation between eignevalues, whose universality has been very 
much discussed over the recent years. Note the dependence in {Np{X))^ of this 
function. This K=l result (P^ ) is indeed equal to (27re"^)fC(Ai, A2), where the kernel 

i^(Ai,A2)is 

Trf\ \ \ sin[7rA^p(A)(Ai-A2)] 

A(Ai,A2J = 7T ^^^ . (37) 

7r(Ai - A2) 

(In the next section we return to the normalizations. It will be explained how the 
extra- factor 27re~^ is cancelled by the normalization constant h^^i ). 

We can now specialize this formula to the moments of the distribution of the 
characteristic polynomial, by letting all the A's approach each other, i.e. letting the 
x's vanish. Before we do that, we should point out that the procedure to obtain 
these moments is in fact subtle. In principle we could have set all the A's equal at 



an early stage of the calculation. If we returned for instance to ( pT]) we might have 

det(e~^'^°^'') 

replaced the limit of — — — - by A{zi, ■ ■ ■ , Z2k) (up to a factor), but then the 

A(Ai, ■ ■ ■ , X2K) 

saddle-point method to obtain the large A^-limit becomes quite problematic. Indeed 
the Vandermonde of the z's at the saddle-point vanishes and it is necessary to go far 
beyond the Gaussian integration. However it is now straightforward to obtain this 
moment from (051). We obtain 



exp-(iVirF(A))F2,,(A,---,A) 

fo Ar MAAi^2exp(-iVis:) /--^c/m^ J^ A^{ui,- ■ ■ ,uk) .„, 

= (27riVp(A)) ^l[ exp-ti^u^) . (3^ 

Expanding the Vandermonde determinant into a sum over permutations, we find 

9 



xE(-l)^""^ ^o^ r.^ n ..■ ••- .o^ . ^n ^' (39) 



1 1 

P,Q (2K - Po-go-1)! ' ' ' (2K - Pk-1 - Qi^-i - 1)! 

in which P and Q are permutations of the integers (0, ■ ■ ■ , i^ — 1). Therefore 



= (_l)^{^-i)/2;^! det i - = ;r!n \ ^' xo (40) 

^ ' o<i,j<K-i {2K - 1 - J - l)\ Y {K + l)\' ^ ' 

and thus finally 

e^V-{NKV{\))F,K{\---A) = (2vriVp(A))^^e-^^ n'^^^^. (41) 

4 Normalizations and Universality 

We have studied in the previous section a Gaussian ensemble of random matrices 
and found that the result ( p] ) for the moment involved {2ttN p{\))^ times a number 
and one would like to see how general is this result, as far as the dependence in the 
density of states is concerned as well as for the normalization. We shall see that this 
behaviour is quite general, and given a proper normalization, that the prefactor is also 
universal. Indeed let us recall how the K-point correlation function of the eigenvalues 
are defined in an ensemble of hermitian N x N matrices X with a probability weight 
proportional to exp — A^TrV"(X). In ^ one finds 

X A2(Ai,--- A^). (42) 

Comparing with our initial definitions (|^) we see that one has the relation 

Rk{Xu---^k) = ^^_-^^, |^^ |exp-iV^V-(A,)|A^(Ai,---A^) 

X F2k{Xi, Xi,- ■ ■ , Xk, Xk)', (43) 



10 



the r.h.s. reduces, up to a normalization, to our previous product of characteristic 
functions of matrices (A^ — K) x (A^ — K), each one beeing repeated twice. On the 
other hand it is well known (|l[]) that this K-point function may be expressed in terms 
of a kernel K^ as 

i?i^(Ai,---,A^) = det irjv(A„Aj), (44) 

and without entering into the precise definition of Kn in terms of orthogonal polyno- 
mials, one should simply recall that —-— ' ,, is universal in the Dyson limit (01) 

p((A + /i)/2)) ^^^ 

(A — /i goes to zero, N goes to infinity, N{X — fj,) finite), i.e. it is independent of the 
polynomial V which defines the probability measure. 

Therefore we define a modified weight, and modified moments. 



^2i^(Ai,A2,---,A2i^) = .^_"^^, ^"'^ sexp— ^^V(A;) \F2k{\iA2,---A2k) 

(45) 
and 

^'^^^^ = {N^~Ky^l7 ^"""^ -^^^(^)> ^2i.(A, A, ■ ■ ■ , A). (46) 

The universality of level correlations implies the universality of M2k- Therefore 
we have to return to the Gaussian case, in order to take into account this new 
normalization, and then the result will be universal. 
lYiom. (|T3|) we have 

N\ Zn_k 1 



{N-K)\ Z^ Ym'-Kh^ 



(47) 



and, given the explicit expression (|T^) of /i„ for the Gaussian case, we find, in the 
large N limit, 

^- ^"-" =(27r)-^e^^. (48) 



{N-K)\ Zn 
With this normalization the universal moment M2k{X) is given by 

M2k{\) = {27r)-^'{2nNp{X)f' YL j^^_ (49) 

11 



In fact this connection between the usual correlation functions and the expec- 
tation values of a product of characteristic functions, (|43|) and (0), allows one to 
recover directly the moment M2k{^), by using the universal expression for the kernel 
K{Xi, Xj) in the Dyson limit, 

/f(A..A,) = -"'':y\-;^)' . (50) 

TT{Xi - Xj) 

The integral representation, over 2K variables describing contours around the K 
poles A/, 

deti<i<j<KK{Xi,Xj)_ 1 J-^ dui r -£- dvi A{ui, ■ ■ ■ , uk)A{vi, ■■■,vk) -^ 



A2(Ai, ■ ■ ■ , Ak) i^! / Y 2m J Y 27rz n.ti Uf=^{u^ - X,){v, - A,) t\ 

(51) 

allows one to write easily the limit in which all the A's are equal: 

deti<i<j<KK{Xi,Xj) 1 fJ^duifJidviA{ui,---,UK)A{vi,---,VK) 



A2(Ai,---,Ax) ir!7V27rz/V2« nL(w. - A)^K - A)^ 

K 

X n^(«^'^0 (52) 

Since the kernel is a Toeplitz matrix, i.e. K{Xi, Xj) = K{Xi — Xj), one can shift the 
m's and the f 's of A and the r.h.s. becomes independent of A. In the case of the sine 
kernel we obtain, in the limit in which all the A's are equal, 

1 I-TT^^i r -1^ dvi A{vi,- ■■ ,vk)A{xi,--- ,xk) -1^ sm{7!-Npxi) 



K\J\^2mJ\^2m a=i[(wi + a;i)^^;f ] r=l ^rx^ 

- (2vr)^ 1}AI + K)V ^''^ 

We have indeed recovered, for any function V defining the probability distribution, 
the universal moment fE9l 



5 Large N asymptotics 

Rather than starting, as in the previous sections, of exact expression for the correla- 
tion functions of characteristic functions, and at the end letting A^ go to infinity, we 

12 



may use a different method to investigate directly the large N limit for the moments 
of their distribution. This method applies for a general probability distribution of 
the form @ and it may also be used to the more general case of an external matrix 
source coupled to the matrix X |^ in this distribution. It turns out that here again 
it is neccessary to consider first F2K for different A/s, and let go all the A^'s approach 
the same A at the end of the calculation. 
/^From (^, we have 

^^ = MG.(A.) (54) 

where Gx{\i) is the resolvent, 

The bracket here denotes an expectation value with a weight which includes both 

2K 

P{X) and Y[ det{\i—X). We assume that the asymptotic spectrum of the eigenvalues 

1 
Xi of X fill a single interval [a, (3] in the large M limit. (It is sufficient to consider the 

single cut case, since we are interested in Dyson short distance universality, which 

involves only the local statistics). Therefore Gx{z) is also analytic in a plane cut from 

the interval [a,/5], and 

Gx{x±tO)=Gxix)Tircpx{x) (56) 

where Gx{x) = [Gx{x + iO) + Gx{x — i0)]/2. The saddle point equation in the large 
M limit becomes 

2K 1 

2MGx{z) - NV\z) + Y. = 0. (57) 

The last term of (0) is of relative order 1/A^ and thus we have to solve this Riemann- 
Hilbert problem to this order. At leading order, we have 2G{x) = V'{x), and up to 
order 1/A^, 

1 2K 

Gx{z) = G{z) + j,{Ga{z) + E^a.(^)). (58) 

lYioia the saddle point equation (|57|), we have Cg{x) = {N — M)G{x) and 

Cx.(x) = -— -. We now set M = iV - K. The functions Gg{z) and Gx(z) 

2(Ai -x) 

13 



are uniquely determined from their analyticity in a plane cut from a to /9, and their 
fall-off as \j z^ for large z (since both G\(^z^ and G(^z^ behave as Xj z at infinity). The 
result is 

Cg{z) = KG{z)- ^ 

^{z-a){z-f3) 

1 1 .^ ^{z-a){z-(3)~^{\-a){\~(3) 



CxAz) 



^iz-a)iz-(3) 



z — A,- 



(59) 



These expressions lead to 



{N-K)Gx{X^) = NG{X,)-^^\og^{K-a){X,-f3) 



1 2K 



y- - Q JXj - a){\j - f3) 

2^.i-^,A.-A/ ^(A,-a)(A.-/3)^- 



(60) 



Since there is a branch cut between a and f3, one must specify whether Aj approaches 
the real axis from above or from below. The sign of the square root on both sides of 
the cut will be denoted e^. There are then a priori 2'^^ saddle points corresponding 
to the different choices of e,-. For each choice of the e,'s, we have 



-J ^1 



r\ 1 ^ 

— logF, = eiMnpiXi) + --— log ^{\i -«)(/?- A,) 



If 1 .^ 6,V(A,-a)(/J-A,) 

where F^ means the value of F2K for given e/s multiplied a factor exp(— ^ Z) V{Xi)). 
Introducing the parametrization 0(a;), defined by a; = ^{a + (3) — ^{(3 — a) cos(f){x) 
and |(/3 — a) sin0(a;) = J{x — a){f3 — x), we have 

d , ._^ei(f){Xi) - ej(j){Xj) 



dA-^°S^^"^ 2 



,V'(A. - a){P - A,) + 6,7(A, -a)iP- A,) 



y^(A,- «)(/?- A,) A, - A, 

(62) 
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Thus we obtain F^ by integration, 



2K „• / ei'/'(Ai)-£jt/'(Aj) N 2K 



^ .A, 

II exp(ejiA^7r / p{x)dx) (63) 

.■_i ^An 



X 



We have to sum over all the saddle-point contributions, i.e. sum over all the different 
choices of ej's. We focus now on the Dyson limit in which the differences Aj — Xj are 
all of order 1/N. Among the 2'^^ possibilities, we retain only the j^ solutions 
in which half of K among e; are positive, and the remaining halves are negative. 
Otherwise, the exponential factor in the final result gives very rapid oscillations in 



the large N limit. This situation is thus exactly similar to that of the previous section 

Agai 
integral 



Again the sum over the ^ saddle-points is conveniently written as a contour 



When we set all the Xj = X, this becomes 

^ = ^/ (2^)- Ug^v^- cos(2|:.„vrArp) (65) 



and we recover the result (|38|). However in this method, since we re-integrated the 
logarithmic derivative of F2K, the constant of integration remains undetermined. We 
may fix this constant by the same requirement that we have used in the previous 
section, and the final result agrees then with the previous calculation. 

6 Symplectic group Sp(N) 

We have studied up to now unitary invariant measures, characterized for the proba- 
bility law of the eigenvalues by the factor |A(a;i, ■ ■ ■ ,xm)|^- We could also consider 
the Gaussian orthogonal ensemble (GOE, with /? = 1) or Gaussian symplectic (GSE, 
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with /? = 4). If we took the GOE for instance, we could immediately relate the corre- 
lation functions of characteristic determinants, to the correlations of the eigenvalues, 
as in (|43| ) (except that since P is one no doubling of the A's is needed), and therefore 
relate the moments universality to the Dyson universal limit. Remaining still with 
the unitary (3 = 2 class, in Cartan's classification of symmetric spaces, we find ensem- 
bles which are invariant under Sp{N) or 0{N). One of the physical applications of 
random Sp{N) matrices, is the statistics of the energy levels inside a superconductor 
vortex P]. In number theory, it is known that some generalizations of Riemann's 
C-functions, such as Dirichlet L-function L{s, Xd) where Xd is a quadratic Dirichlet 
character of mod \d\, present a spectrum of low lying zeros on the line Re s = 1/2, 
which agrees with the statistics of the eigenvalues of the Sp(N) random matrix theory 
121 , |l^. In this Sp(N) invariant symmetric spaces, the eigenvalues appear alway in 



pairs of positive and negative real numbers. Due to this fact, a new universality class 
governs the correlations of the eigenvalues near the origin, i.e. near s = 1/2, (whereas 
in the bulk one recovers the previous unitary class). 

Therefore we study now the new universality class, which governs the new scaling 
near the origin. We thus consider random Hermitian matrices X, which are 2M x 2M 
and satisfy the condition 

X^J +JX = (66) 

where J is 

The unitary symplectic group is a the subgroup of SU{2M) consisting of 2M x 2M 
unitary matrices, satisfying the symplectic constraint 

U^ = -JUU (68) 

The integration over this unitary symplectic group for Fk{Xi, ■ ■ ■ , A^) gives 

K 



a=l 
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1 „ M M K M 

Repeating the analysis of section 2, Fk{Xi, ■ ■ ■ , A;^) is given again by a determinantal 
form as (]I^. Changing xi to x^ = yi and denoting /ij = A^, we have 

.oo K K ^ K K 

F^(/ii, ■■■,^iK) = J UdmU yl I[iy^ - Vif IT IK/^- - i/^e"^^^^ (70) 

1=1 i=l i<j a=l i=l 



The orthogonal monic polynomials for this measure are the Laguerre polynomials 

(-) 
Ln (y), which is defined by 



;-i)",,, /(i!<(i + 'i.)"+4,_„„. 



(-) 
normalized as required to Ln^ (y) = y"' + lowerdegree. The orthogonality condition 

is 

j^" dye-''^^d^\y)Lll\y) = hJn,m (72) 

with hn = n\T{n + |)/A^^"+2, and Hn-i — 27re"^^ in the large N limit. 
From (|I^) , we have 



j^KiM+f-^) A(/i) 
K ^-y ^ (^ ^ ■y\M+K-^ K -y r 

i=l ^"^ /=1 ^; j^j i -t- ^i J- -t- ^J 

X-l 

We now set M = AT - ir, and the factor [] (^ + /)!/Ar^(*^+-^/2-i/2) -^ ^^^^j ^^ 



{2TiN)^e~^^ , up to corrections of relative order 1/A^ in the large N limit. The large 
N limit is governed by the saddle-point equations zf + zi + — = 0. In the following 
we study the scaling vicinity of the origin, in which all the /i;'s scale as l/N"^. Then 
zf at the saddle-point may be expanded 

zi^^-l + 0{^,) (74) 
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where e; = ±1. 

Noting that m-r^ - T^) = IK-^'-"^ + ^>j) ^ ^''^'''''>^'I[i^^^^ ' ^3^) 
(cj = ±1), and combining it with the Vandermonde A(A^), we are left with a fac- 
tor Yli<:j ^_x +e X ^^ ^^^^ scahng hmit. We have also the exponential e~^^^"^" = 

We have again to sum over all the saddle-points, which are characterized by 
the sign of ej = ±1, and to include the factor due to the fluctuations near the 
saddle-point. The Gaussian fluctuations yield a factor (2n/(—2iei\^N)y^'^. Then 
{1 + Zi)"^^"^ ~ (Aj/(eji))^/^. There is a l/(27ci)^ in addition. We have an extra Cj due 
to the contour direction, which goes through two saddle points ; one is in the positive 
imaginary plane and the other in the negative half-plane. When K=2, and Ai and 
A2 are of order 1/A^, we obtain 

27re-2^ 
i^2(Ai, A2) = -——Ksp{\u A2), (75) 

A1A2 

with the kernel Ksp{\i, A2) given by, 

sin[iV(Ai - A2)] sin[iV(Ai + A2)] ._. 

^"^^^^' ^'^ - 2.{X, - A2) 2.{X, + A2) ^^^^ 

The coefficient (27r)e^^^ is cancelled by the normalization factor 1/hN-i- Putting 

Ai = A2 = 0, we have neglecting the factor 27re^^^, ^2(0) ~ ^^^^• 

For general K, Fk{Xi, ■ ■ ■ , A^) becomes in the scaling limit 

Fk(Ai,---,A^) = (-l)^(^-^+^)(2vriV)fe-^^«f(^-^)(^)f ^ 



The sum over all the saddle-points, characterized by e^ ± 1, is conveniently written 
as a contour integral. 



y 



l)f (^-)^ / . . . / n(— ) .^y^^"^ e-^gt."') (7^ 
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where the contour encloses Ui = ±Aj. We may now set Xj = A, and keeping track 
of various coefficients, we obtain the K-th moment Fx{X, ■ ■ ■ , A). For general A, the 
result has a complicated form, but when A = 0, it becomes a number 



r)k/2^-NK 

Fk{0,---,0) = ^4 ivf(^+i)(z)f(^-3)(-l)^(^-i) 

k\ 

This representation allows one to compute the K-th moment at the origin. By the 



Expansion of the VanderMonde determinants, similarly to (|39D, (|79D is reduced to a 
determinant form. We have by the normalization; Fk{0) = (27r)^~e^^FR-(0), 

FK{0) = {-lf^U7^,^-^ (80) 

Comparing to the result of the unitary case in (^), we notice that the exponent of 
N is different and the universal coefficient is given also by the product of the ratio of 
the factorizations. 

For F2k{Xi, Ai, ■ ■ ■ , Xk, Xk), the even 2K-th moment may be obtained again from 
F2k{Xi, Ai, ■ ■ ■ , Xk, Xk) = det[ii'(Ai, Aj)]/(A^(A^) Jl Af) ; using the expression for the 
kernel (0), we have for the 2K-th moment, 

det[K{Xi,Xj)] 2^ r ^ dui r ^ dvi A{u^)A{v^) 



i^^nem 



A2(A2) n A? K\J^^2mJ^^ 2m Ul, nf=i(nf - A,^) Ut^ UUi^f - X]) 

1 -^sinlN (ui-Vi)] ,^ 



(27r)^ f=i Ui-Vi 

For general A, the result has a complicated form, but again here one can compute 
form there the values at A = 0. The result agrees with the previous expression of 
-^2i^(0) in (^Ol). One may also use the large N asymptotic analysis as in section 5 and 



rederive the results as the same sum ([78|) over the saddle points. 
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7 Orthogonal group 0(N) 

We discuss here the 0{2N) case, which is different from Sp{N) ( whereas 0{2N + 1) 
has a structure which is similar to Sp{N) 0). In number theory, for example the 
twisted L-function, L^(s,Xd) presents a spectrum of low lying zeros, which agrees 
with the statistics of the eigenvalues of the 0{2N) random matrix theory |jl2|, |14 . 
Then in terms of eigenvalues 

K 

Fk{Xi,---Ak) = <ndet(A«-X)> 

1 „ M KM 

= -^ ]1M^^) ^\4r--,xl,)l[l[{Xl-x',) (82) 

The difference between the symplectic and orthogonal case is due to the absence 
of the factor l\x^. Using the analysis of section 2, Fk{\i, ■ ■ ■ , Xr) is given by the 
determinantal form as in (|I1|). Changing x] to i/i and denoting ^i = \f, we have 



l-oo K K -^ K K 



(--) 
The orthogonal polynomials for this case are Laguerre polynomials Ln ^ (y), which 



is defined by 



Li.-''(!/) = (-l)"^f^(A)»(y..4e-«>') 



'K, I ^dl^e-"" (84) 



{--) 
normalized as Ln ^ (y) = y^ + lowerdegree. The orthogonality condition is 



r dye-^'y^Li '^\y)LL~'\y) = hj„„m (85) 

Jo ./y 



Vy 

with hn = n\T{n + |)/A^^"+2, and h^^i — 27re^^^ in the large N limit. From (jT^) 
we have similar to the Sp{N) case. 



/M^2.r.\y .r^ ^4^4+.. 1+., 






(86) 
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K-1 

KN 



We set M = N-K, and the factor [] (M+Z)!/A^^(^+^/2-^/2) is equal to (27rA^)^e 



, up to corrections of relative order 1/A^ in the large N limit. The saddle point zi is 
same as (JT^)- The only difference is the extra factor (1 + zi)^^ ~ j^. When K=2, we 
obtain 

F2(Ai, A2) = 27re-2^iro(Ai, A2), (87) 

with the kernel Ko{Xi, X2) given by, 

sin[iV(Ai - A2)] , sin[iV(Ai + A^)] 
"^^^^^'^^^^ 2.(A,-A2) + 2vr(A, + A2) ^''^ 

The factor (27r)e^^^ is cancelled by the normalization factor l/hj^^i ~ (27r)^^e^^. 
Putting Ai = A2 = 0, we have neglecting the factor 27re"^^, ^2(0) ^ ■^{2N). 
For general K, Fk{Xi, ■ ■ ■ , A^) becomes in the scaling limit 

FKi\i,---AK) = (-l)^(^-^+^)(27riV)fe-^^(z)T(^-3)(^)f 1 



X Yl V^ TT (89) 

The sum over all the saddle-points, characterized by e^ ± 1, is conveniently written 
as a contour integral, 

- <-'*''--4/--in(S)|^ff£f e-EL...) ,00, 

where the contour encloses Ui = ±Aj. We may now set \j = A, and keeping track 
of various coefficients, we obtain the K-th moment Fx{X, ■ ■ ■ , A). For general A, the 
result has a complicated form, but when A = 0, it becomes a number 

nk/2^~NK 

F^(0,...,0) = ^-4 ivf(^-i)(z)f(^-5)(_i)^(^-i) 

k] 
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The normalization factor is (27r) 2 e^^ for Fk{X). Denoting the normahzed K-th 
moment by Fk{X), we have 



K-l 



n (2iv)f(^-i) 



Fk{0) = (-1)^^ n T^^T (92) 

We have F2k{^i, -^1, ■ ■ ■ , Xk, Xk) = det[K{\i, Aj)]/A^(A^). Using the expression for 
the kernel (0), we obtain for the 2K-th moment in the orthogonal 0{2N) case, 



A2(A2) K\T^^2mT^^ 2m Ul, nf=i(«? - A,^) nf=i U^i^f - \]) 

^^ 1 ^ sin[iV(n,-t;,)] ^ ^^^^ 



(27r)^ j=i Ui - V: 
Inserting Aj = 0, we find the consistent result with (pl|). 

8 Negative moments 

In the number theory literature one finds various moments in which powers of the 



zeta-functions appear in the denominator [|Tl|]. The equivalent for random matrices 

K 

would be to consider expectations values of the form < JJdet(A/ — XY'- > in which 

1 
the e's are ±1. One cannot use any more the techniques introduced hereabove but, 

at least in the Gaussian case, it is easy to obtain exact expressions through the use 

of auxiliary integrations, over both commuting and anti-commuting variables. 

We first rederive our previous results for positive moments (i.e. ei = +1 for all 

Vs) . Let us introduce M Grassmann variables c^, c^ and an integration normalized 

to 

r dcdc , ,, 

/ cc = 1. 94 

J vr 

Then, for an hermitian M x M matrix X, one has 

r *^ dr dr 
det(A-X) = iV-*^ /n^-^exp^iV^[c-(A5,,,-X,,,)c,]. (95) 
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K 



A product Y\^ det{Xi —X) is represented by a product of K integrals of the type (|9 

1 
At the end the random matrix X appears in an expression of the form 

K M 



exp-^AT^ 5: X.,c«c?. (96) 



1=1 a,b=l 



With the Gaussian probabihty weight (0) we have 

< exp iA^TrAX >= exp -— TrA^ (97) 

and thus 

< n det(A, - X) >= N~'' / n ^^^^^ exp (AT Y. tXau + ^ E 7/,n7w) (98) 
with 

N 

llm = Y.'^^I^C^:'^. (99) 

a=l 

We can use an auxihary K x K hermitian matrix B to replace the quadratic terms 
in 7 by 

exp — Tr72 = (—)KV2 f ^k^ ^ ^^^ {NTi-fB - -TrB^). (100) 

2 2rr J 2 

We are left with an integral over the Grassmannian variables 



■' a=ll=l '''^ l,m=l a=l 

det {\i6i,^-tBi^)) . (101) 

l<l,m<K / 

We end up with an integral over a. K x K hermitian matrix B : 
< n det(A, - X) >= i^f"^' J d^'B {det{XAm - tBi^)}'' exp -yTrB^. (102) 

Therefore, from this method as well, we have reduced the correlations of the charac- 
teristic functions of the matrix, to an integral over K^ variables. If one is interested 
in the moments, i.e. A/ = A for all /'s, one may take as variables the eigenvalues 6/ 
of B (which yields a factor A^(Ai, ■ ■ ■ , A;^ )), and recover the previous expressions. 
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For the A;'s non-equal, one must first shift the matrix B of the diagonal matrix 
i^i^im, and then integrate out the unitary group SU{K) by the Itzykson-Zuber for- 



mula |jT5|, [T^, 0, to reduce it, as before, to an integral over K variables (a slightly 
different integral, but which may be handled in the large iV-limit in an identical 
fashion) . 

In case of negative moments the method is identical, except that we need now 
ordinary commuting variables, instead of Grassmannian. Indeed starting from 

__J__ ^ ^M / n ^ exp ±m urns., - X., ± ^.S^M (103) 

one can introduce, for each factor (det(A/ — X)y' an integration over M complex 
variables {(j)l,(j)a) if ^i = —1, or over M complex Grassmannian variables {ca,Ca) if 
ei = +1. The expectation value with the Gaussian weight P{X) is then immediate. 
Of course for the negative moments, one must pay attention to the sign of the in- 
finitesimal imaginary part of the A's since there is a cut on the real axis along the 
support of Wigner's semi- circle. 

Although the method is obvious and elementary, the notations can become cum- 
bersome and, rather than working out the most general case, and arbitrary choices 
for the signs of the imaginary parts, we restrict ourselves to an example. If we con- 
sider only negative powers, we may follow identical steps as hereabove with positive 
powers, and we find 

1 (let I A; — A -\- 16) Z7T J 



X exp TrBl (104) 



When all the A's are equal the r.h.s. simplifies to an integral over K variables 

/ W^hi — jT-. — '- ^ — -ryexp— ( — y^&?)- For the A^'s non equal, after a shift of 

J 1 111 (A - 6i + leY^ 2 ^ 

the matrix B and the integration over SU{K), one obtains 

^\idet(A,-X + ^6) ^"^Tr^ '^''P 2^^' 
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x'n „i,.fg::::::S exp-^i:(i.?H-.,A,). im 

from which one could repeat easily the analysis of section 3. 

9 Discussion 

We have discussed the universal expressions for the moments of the characteristic 
polynomials in a random matrix theory, where the ensembles belong to the unitary 
family {(3 = 2). 

We have shown that these universalities are related to the universality of the kernel 
in the Dyson's short distance limit. Since the statistics of the zeros of the (^-function 
follows the universal behavior of Gaussian unitary ensemble (GUE) [|12], |T^, the power 
moment of the (^-function also has to follow the universal behavior of GUE. We have 
studied here the characteristic polynomial, which corresponds to the ^-function on 
the critical line, and we have found a universal behavior for the moments of the 
characteristic polynomial. The universal number (|49|) appears indeed in the average 
of the moment of the (^-function, which was conjectured as (H), 'Jk = Uo^^ l\/{K+l)\. 

Our method of the splitting the singularity by the introduction of the distinct 
Aj may be applied directly to the average of the power moment of the Riemann 
C-function. We consider the average of the product of (i^i)^ Si = ^ ± i{Xi + t), 

1 r^'^^ 

p=T nc(^o^^ (106) 

where we choose K positive Aj's and K negative ones. If, at the end of the calculation, 
we set all the positive Aj's equal to A and the negative ones to —A, one recovers the 
2i^-th moment of the modulus of the (^-function. When T is large, the leading and 
the next leading terms of the derivative of InF with respect to Aj, are presumably 
given by 

— -~±nnT-5:- -, (107) 

OAi ^y • Ai - Aj 
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where the pole in the second term appears when two distinct Aj coincide, one of 
the Aj's with a plus sign and the other one with a minus sign. In the Appendix a 
discussion of the assumptions leading to (|107| ) is given. Then, after integration, we 
have, following a line of arguments similar to those of section 5, 

Therefore, if we let the Aj's coincide, we recover the integral (^), which provides the 
universal coefficient ''Jk- The coefficient c is not determined by this method, which 
starts with the logarithmic derivative of F, and an extra normalization condition 
is needed. In (A. 7) it will be argued that a coefficient qk is present in the result, 
which is the residue at s = 1 of a function gxis) defined in the Appendix ; it is thus 
plausible that the coefficient c in ( |108| ) is nothing but c = ax- 

We have also investigated negative moments as ( |105| ). This result may apply to 
the mean value of negative moments of the (^-function. Indeed, the exponent K^ of 
logT for the negative integer K, has been conjectured ||TT[] . 

For the symplectic and orthogonal case, Sp{N) and 0{2N) ensembles, there may 
be also be a correspondence between the random matrix results (0), (0) and the 
average values of the certain L-functions, with the same jk, as far as there is a uni- 
versality. Existing conjectures [|I^] for the moment of the L function shows the same 
exponent y(^ + ^) ^^^ t(-^ ~ ^) ^'^^ ^^^ symplectic and the orthogonal cases, how- 
ever, the conjectured values of jk is different from our result (^) for the symplectic 
case. 
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Appendix : Summation formula for the Riemann zeta-function 
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The Riemann (^-function is given by 

n=l "" P ^ p^ 

where p is a prime number. The i^-th power of this function is written as 



oo 



[C(^)]^ = E 



n=l n 



p p p 

where dxin) is the i^-th Dirichlet coefficient. When ra is a power of the prime number, 
dxip') = ^{K + j)/r{K)j\, (this follows easily from the definition of the Dirichlet 
coefficient dx{n) = E -'-)• 

ni---nx=" 

We consider now the average of ( |A.2|) on the critical line s = ^ + it over a large 
interval T, 

1 r^ ..,..1 ...,o^ 1 r^ ,..^dK(n 



I dt\a^+u)r = u dt\j:^\^. (a.3) 

Jo 2 T Jo ^1 722+** 



T Jo Z 1 JV „^^ 77,2 

oo ^ / \ 

Expanding the sum I T^ P, which appears in (|A.3|), we first examine the diag- 

n=l "- 

onal terms , 

= na-p'r^^a- '^'^'^4~'V " + ---) 

= [as)f'9K{s), (A.4) 



where 



.2 °^ -^2 



^i.(^) = nia-p"r^E^]- (A.5) 

P j=0 P 

The function gxis) is an analytic function of s, including the point s = 1. 

Let us examine the contribution of these diagonal terms given by ( [A.4| ) to ( |A.ci| ). 
Their contribution is conveniently found, if we apply the following inversion formula 
(Perron formula). 

oo 
B{S) = Y.^nU-' 
n=\ 
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fix) = E^n (A-6) 



n<x 

Then, we have 

1 rc+ioo 

fix) = — / Bis)x's''ds, (A.7) 

in which c is some arbitrary real positive number. Substituting 6„ = d\in), and 
Bis) = C^ is)gKis), we obtain, from the residue of the singularity at s = 1, 

E4H = ^x\og^'-'x + Oix\og^'-'3x). (A.8) 

By a partial summation, this approximate calculation yields , 

where qk = 5'x(l)- 

From these formulae, it is seen that the contribution of the diagonal terms to 
the average of the K-th power moment of the (^-function does take the asymptotic 
form of (^. However, neglecting the off-diagonal terms, we failed to reproduce the 
proper coefficient 7x, whose understanding clearly requires the off-diagonal products 
in ( |A.4|) as well. 



A lower bound for the 2i^-th moment is known 18 



rT+Y 1 „ 

/ \Ci- + it)\^^ dt ^ Y \og^ Y (A.IO) 

Jt-y 2 

where log*^ T < Y < T. An upper bound seems difficult to obtain, and (0) remains 
as a conjecture, except for the K = 1 and K = 2 cases,for which it has been derived. 



We note here the results and the conjecture of Montgomery [jT9l about the density 
of the zeros of Riemann ^-function and their correlation. When 7 is a zero on the 
critical line, (i^ + i'j) = 0, 

E l>(| + o(l))f logT (A.ll) 

0<7<T "^ ^^ 



pj3 sin TT?/ 

E l = (l + o(l))[/ (l-( )')du + 6ia,l3)]TL (A.12) 
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where L = logT / {2tt) , and 6{a, /?) = 1 for G [a, P], and otherwise zero. Then ( |A.11| ) 
is equivalent to the average density of state in (|ID, with for K = 1, ■jk = ciR = ^ and 
(|A.12|) is equivalent to the pair correlation function in random matrix theory. 

Les us present the arguments which lead to the conjectured formula (|107| ) ; we first 



assume assume that Ai — A2 ~ 0((lnT)~^) for large T. The diagonal approximation 
for the product of ({si) and ({^2), which earlier gave the expected behaviour for 
the moment, but with a wrong coefficient, may thus be applied here again, since we 
are taking a logarithmic derivative, which is unsensitive to overall normalizations. 
Within this assumption, we obtain 



n<T 



dXi 



ln[ / dx- 



dXi 7to x1+^(^i-^2)J 

~ -ilnT-- (A.13) 

Ai — A2 

We have considered up to now what happens when Ai — A2 is small, but we should 
repeat the same arguments for the Dyson limit in which all pairs Ai — Xj are of order 
(logT)"^. Therefore, when one sums over all possible combinations, one obtains 
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